ZU-TH 15/06 



1 



A numerical approach to the double real radiation part of e + e — ► 3 jets 
at NNLO 

G. Heinrich a * 

a Institut fiir Theoretische Physik, Universitat Ziirich, Winterthurerstrasse 190, 8057 Zurich, Switzerland 
We report on the sector decomposition approach to the double real emission part of e + e~ — » 3jets at NNLO. 



1. INTRODUCTION 

Precision measurements at high energy colli- 
ders in the recent past led to stringent tests of the 
Standard Model and important bounds on New 
Physics. Measurements of jet rates and shape ob- 
servables in e + e _ annihilation are of particular 
importance, as they allow for a precise determi- 
nation of the strong coupling constant a s . At 
hadron colliders, a cross section involving n jets 
is proportional to a™ at leading order, such that 
an accurate knowledge of a s will be important at 
the LHC, where many interesting processes con- 
tain jets in the final state. 

A determination of a s from jet rates and shape 
observables is described e.g. in pQ, where one 
can see that for LEP measurements, the experi- 
mental uncertainty is smaller than the theoretical 
one, which is based on resummed next-to-leading 
order (NLO) calculations. As the theoretical er- 
ror is dominated by scale uncertainties, an NNLO 
calculation will considerably improve this situa- 
tion. A future International Linear Collider will 
allow for precision measurements at the per-mille 
level, which offer the possibility of a determina- 
tion of a s with unprecedented accuracy, provided 
that theoretical predictions at NNLO are avail- 
able. 

2. COMPUTATIONAL METHODS 

The calculation of e + e~ — * 3jets at order a? 
requires the calculation of virtual two-loop correc- 
tions combined with a 1 — -» 3 parton phase space, 
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one-loop corrections combined with a 1 — > 4 par- 
ton phase space where one parton can become 
soft and/or collinear ("unresolved"), and the tree 
level matrix element squared for 1 — > 5 partons 
where up to two partons can become unresolved. 
The unresolved particles lead to a complicated 
infrared singularity structure which manifests it- 
self in 1/e poles upon phase space integration. 
These singularities have to be subtracted and 
cancelled with the ones from the virtual contri- 
butions before a Monte Carlo program can be 
constructed. To achieve this task, two differ- 
ent approaches have been followed, one relying 
on the manual construction of an analytical sub- 
traction scheme [2l3l4l5lfil7l8| . the other one re- 
lying on sector decomposition [91101111121131141 
Il5llrill7ll8ll9l20l2l| . The main features of the 
methods based on the explicit construction of a 
subtraction scheme are the following: The sub- 
traction terms are integrated analytically over the 
unresolved phase space, such that the pole coef- 
ficients are obtained in analytic form. This re- 
quires appropriate phase space factorisation and 
subtraction terms which are simple enough to be 
integrated analytically in D = 4 — 2e dimensions. 
This method naturally leads to a close to minimal 
number of subtraction terms, and allows insights 
into the infrared structure of QCD. 

In the sector decomposition approach, the poles 
are isolated by an automated algebraic procedure 
acting in parameter space, and the pole coeffi- 
cients are integrated numerically. The advantages 
of this approach reside in the fact that the extrac- 
tion of the infrared poles is done by the computer, 
and that the subtraction terms can be very com- 
plicated as they are integrated only numerically. 
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On the other hand, the algorithm which isolates 
the poles increases the number of original func- 
tions, thus producing rather large expressions. 

The application of sector decomposition to real 
radiation at NNLO first has been presented in 
|12I13| . and the combination of the sector decom- 
position approach with a measurement function 
first has been proposed in |14| . A number of 
NNLO results based on this method have been 
obtained meanwhile |lfll4llfilltiil7llHll*)fcW| . Its 
application to the double real radiation part of 
e + e~ — > 3 jets at order [2] is particularly chal- 
lenging due to the high number of massless parti- 
cles in the final state, which leads to a very com- 
plicated infrared singularity structure. 

3. SECTOR DECOMPOSITION 

The wide range of applicability of sector de- 
composition goes back to the fact that it acts in 
parameter space by a simple mechanism. The pa- 
rameters can be Feynman parameters in the case 
of multi-loop integrals, or phase space integration 
variables, or a combination of both. In the follow- 
ing, the working mechanism of sector decomposi- 
tion will be outlined only briefly for the example 
of phase space integrals, details can be found in 
jllll5ll§l21| . 

The phase space integral of a matrix element 
squared, combined with some measurement func- 
tion J defining a physical observable, typically 
contains "overlapping" structures like 
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where only the dependence on two of the inte- 
gration variables is shown for pedagogical sim- 
plicity. In order to extract the poles in 1/e, the 
singularities for x,y — > need to be factorised. 
Sector decomposition is a way to achieve the fac- 
torisation of this type of entangled singularities 
in an algorithmic way: First the integration re- 
gion is split into sectors where the variables x 
and y are ordered by multiplying with unity in the 
form [Q(x — y) + Q(y — a;)]. Then the integration 



domain is remapped to the unit cube: after the 
substitutions y = x t in sector (a) and x = y t in 
sector (b), one has 
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where the singularities are now factorised. For 
more complicated functions, several iterations of 
this procedure may be necessary, but it is eas- 
ily implemented into an automated subroutine. 
Once all singularities are factored out, they can 
be subtracted using identities like 
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where we recognise the form of plus distributions. 
The result can subsequently be expanded in e, 
such that a Laurent series in e is obtained, where 
the pole coefficients are sums of finite parameter 
integrals which can be evaluated numerically. 

For the numerical evaluation it has to be as- 
sured that no integrable singularities are crossed 
which spoil the numerical convergence. In the 
double real radiation part of e + e~ — > 3 jets, sin- 
gularities which are located in the interior of the 
integration region do indeed occur. However, it 
is always possible to remap them to endpoint sin- 
gularites by a convenient variable transformation. 
After such a remapping, they are amenable to sec- 
tor decomposition. 
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In order to calculate the double real radia- 
tion part of e + e~ — > 3jets at NNLO, sector de- 
composition is applied to extract the poles ap- 
pearing in massless 1 — > 5 particle integrals. 
As a simple example, let us consider the 5- 
particle cut of the ladder graph shown in fig. 1. 
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Figure 1. The ladder graph 
Sector decomposition leads to |23 



TV 



3 /a s \3 f 0.16662 



+ 4 [1.4993 - 0.4999 log [ - 
e 2 VM 

+- [5.5959-4.4978 log f ^ 

e 

+ 0.7498 log 2 + finite! , (3) 

where the numerical accuracy is better than 1%. 
The correctness of the result can be checked by 
exploiting the fact that the sum over all cuts of 
a given (UV renormalised) diagram must be in- 
frared finite. This is shown diagrammatically in 
fig. |21 After UV renormalisation, we obtain the 
condition 

Ti^5 + %i T\^a + z 2 Ti_>3 + z 3 Ti^ 2 = finite , (4) 
where Ti_>j denotes the diagram with i cut lines. 
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Figure 2. Cancellation of IR divergences in the 
sum over all cuts of the renormalised graph 

The renormalisation constants Zi (in Feynman 
gauge) are given by |15I21| 
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The expressions in eq. Q for i < 5 combine to 
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We can see that the poles in © are exactly can- 
celled by the 5-parton contribution (J^J within the 
numerical precision. 

4.1. Differential results 

Although the sector decomposition approach 
is considered to be a "numerical method" , as 
the pole coefficients are only calculated nume- 
rically, the isolation of the poles is an algebraic 
procedure, leading to a set of finite functions for 
each pole coefficient as well as for the finite part. 
These finite functions are written to a Fortran 
program and evaluated numerically using Monte 
Carlo techniques. In order to obtain results which 
are differential in a certain physical observable, 
any (infrared safe) measurement function can be 
included at the level of the final Monte Carlo pro- 
gram, which means that the subtractions and ex- 
pansions in e do not have to be redone each time 
a different observable is considered. Further, the 
measurement function does not have to be an an- 
alytic function, but can be a subroutine acting 
on the four-momenta of the final state particles, 
as it is typically the case for a jet clustering rou- 
tine. The reason why such a maximal flexibility is 
possible resides in the fact that the program de- 
scribed in has the architecture of a partonic 
event generator. This requires to keep track of the 
mappings of the original phase space variables to 
the variables used in each sector, i.e. each end- 
point of the iterated decompsition tree. These 
mappings will be different for each sector (cf. the 
arguments of the function J in the first and sec- 
ond line of eq. 10)), but keeping this information 
allows to express the energies and angles of the 
final state particles in terms of the sector vari- 
ables and thus to reconstruct the fully differen- 



4 



G. Heinrich 



tial information on the final state. In this proce- 
dure, function evaluations in sectors which do not 
pass the kinematical requirements imposed by the 
measurement function are unavoidable. Having 
to deal with a large number of sectors, this can 
lead to a serious drop in efficiency of the Monte 
Carlo program. Therefore, in order to construct a 
program which produces results within a reason- 
able time scale, it is crucial to keep the number 
of sectors low. This can be achieved by (a) using 
optimised phase space parametrisations for each 
topology, (b) using information on physical limits 
in the decomposition algorithm. 

Optimising the phase space parametrisations 
means the following: The matrix element squared 
can be divided into a certain number of topolo- 
gies, which are defined by a certain set of denom- 
inators (which will be combinations of Mandel- 
stam invariants, see e.g. eq. In a given 

phase space parametrisation, some of the invari- 
ants will naturally be in a factorised form, but 
such a form cannot be achieved for all invariants 
simultaneously. An optimised parametrisation is 
one where the number of non-factorising invari- 
ants in the denominator is kept minimal. There- 
fore one has to choose different parametrisations 
for each topology, resp. for each class of topolo- 
gies with the same factorisation properties. 

Using information on physical limits exploits 
the fact that in most cases of physical relevance, 
the measurement function is such that it would 
prevent certain poles from arising at all if it were 
included in the e-expansion. For example, poles 
associated with a 2-jet configuration, where 3 of 
the 5 final state particles become theoretically un- 
resolved, will be killed by a 3-jet measurement 
function. It would therefore be desirable to sup- 
press the terms associated with such configura- 
tions already at the level of the e-expansion, in 
order to avoid an unnecessarily large number of 
terms associated with the isolation of these poles. 
On the other hand, we would like to keep the flexi- 
bility to include any measurement function only 
at the stage of the final Monte Carlo program. Fo- 
cusing on the process e + e~ — > 3 jets, this dilemma 
has been solved by including some "preselection 
rules" in the e-expansion which reject configura- 
tions which will surely be 2-jet configurations. In 



the example shown here, this can be achieved by 
introducing a cut parameter y th - which must 
be smaller than any possible experimental reso- 
lution parameter y cut - for the variable S1345, as 
S1345 — > always corresponds to a 2-jet configu- 
ration. In this way, one can reduce the size of 
the expressions considerably without loosing the 
flexibility to specify different jet algorithms later. 

To illustrate the action of the jet function, 3-, 
4- and 5-jet rates using the JADE algorithm 
are shown in Fig. [21 based on the toy matrix ele- 
ment built from the graphs shown in Fig. 




3 jets 

4 jets 

5 jets 
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Figure 3. 3-, 4- and 5-jet rates at order a, for 
the toy matrix element 

As a more complex example, the double real 
radiation part of a non-planar topology (see fig. 4) 

alSO haS been Calculated, magnitude of jet rates tor non-planar topolo 




Figure 4. (a) non-planar topology 

(b) 3-, 4- and 5-jet rates at order al 
for the non-planar topology 



In this case, square-root terms in the denomi- 
nator are unavoidable, which implies that the ex- 
pressions produced by the sector decomposition 
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will be larger due to the presence of more non- 
factorising denominators. The virtual corrections 
to this topology have not yet been included. 

Note that the results shown in figs. |3| and 4 
are unphysical, as they do not contain all contri- 
butions to form a gauge invariant quantity. The 
purpose of these figures is merely to demonstrate 
the action of the jet function and thus the power 
of the method to produce differential results. 

The CPU time is 0(2h) for a precision better 
than 1%. Note that all topologies can be calcu- 
lated in parallel, such that the CPU time for the 
full double real radiation part will be of the or- 
der of the one needed for the most complicated 
topology. 

What remains to be done, besides the inclu- 
sion of the remaining topologies, is the combina- 
tion with the one-loop plus single real and the 
two-loop virtual corrections. How to proceed ef- 
ficiently by combining sector decomposition with 
analytical results is sketched in |17I23| . 

5. SUMMARY AND OUTLOOK 

The method based on sector decomposition to 
calculate e + e~ — * 3 jets at NNLO has the ad- 
vantage that the isolation of the infrared poles is 
done by an automated routine and that it does 
not require the analytical integration of subtrac- 
tion terms. However, the method produces large 
expressions, which is an issue for a process like 
e + e~ — > 3 jets at NNLO where the matrix ele- 
ment to start with is already large. It has been 
described how to limit the proliferation of terms 
in the course of sector decomposition, and differ- 
ential results have been shown for a subpart of 
the process. 

The inclusion of massive particles within this 
method is certainly more straightforward than in 
analytical approaches. In fact, as the masses act 
as infrared regulators, the number of decompo- 
sitions will be less and therefore the produced 
expressions should be smaller than in massless 
cases. 
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